A TRMM rainfall estimation method applicable to land areas 
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Abstract 


Utilizing multi-spectral, dual-polarization Special Sensor Microwave Imager 

(SSM/I) radiometer measurements, we have developed in this study a method to 
retrieve average rain rate, Rf , in a mesoscale grid box of 2° x3° over land. The key 

parameter of this method is the fractional rain area, fR, in that grid box, which is 
determined with the help of a threshold on the 85 GHz scattering depression 

deduced from the SSM/I data. In order to demonstrate the usefulness of this 
method, nine-months of Rf R are retrieved from SSM/I data over three grid boxes in 

the Northeastern United States. These retrievals are then compared with the 
corresponding ground-truth-average rain rate, R g , deduced from 15-minute rain 

gauges. Based on nine months of rain rate retrievals over three grid boxes, we find 
that Rf^can explain about 64 % of the variance contained in Rg. A similar 

evaluation of the grid-box-average rain rates Rgscat and RsrL/ given by the 
NASA/GSCAT and NOAA/SRL rain retrieval algorithms, is performed. This 
evaluation reveals that Rgscat and Rsrl can explain only about 42 % of the 
variance contained in R g . 

In our method, a threshold on the 85 GHz scattering depression is used 
primarily to determine the fractional rain area in a mesoscale grid box. Quantitative 
information pertaining to the 85 GHz scattering depression in the grid box is 
disregarded. In the NASA/GSCAT and NOAA/SRL methods on the other hand, 
this quantitative information is included. Based on the performance of all three 
methods, we infer that the magnitude of the scattering depression is a poor indicator 
of rain rate. Furthermore, from maps based on the observations made by SSM/I on 
land and ocean we find that there is a significant redundancy in the information 



content of the SSM/I multi-spectral observations. This leads us to infer that 
observations of SSM/I at 19 and 37 GHz add only marginal information to that 
given by 85 GHz scattering depression. 

As with other methods, the area-average rain retrieval method developed in 
this study needs tuning with radar and/or rain gauge observations. In the TRMM 
mission, the microwave radiometer rain retrieval algorithm can be tuned with 
TRMM radar observations. Since the radiometer has about 3.5 times wider spatial 
coverage compared to the radar in the TRMM mission, such an algorithm can be 
useful to extend geographically the rain information provided by the TRMM 
Precipitation Radar. 



1. Introduction 


In recent studies, rain rates retrieved from the Advance Microwave 
Precipitation Radiometer (AMPR) on a footprint scale using theoretical models 
(Heymsfield et al., 1996 and McGaughey and Zipser, 1996) showed poor agreement 
with radar rain observations made from aircraft. These radiometer and radar 
observations had a field of view (fov) of a few kilometers. Such poor agreement is 
also observed in the rain rates deduced theoretically from Special Sensor Microwave 
Imager 1 (SSM/I) radiometer observations that have a much larger footprint 
(~ 15 km). During the Global Precipitation Comparison Project (GPCP), SSM/I rain 
retrievals on a footprint scale were deduced from several algorithms and compared 
with ship-borne radar rain observations in the TOGA-COARE 2 experiment. Ebert et 
al. (1996) presented a detailed account of these comparisons. From that GPCP study, 
we infer that rain rate deduced from the algorithms on a footprint scale over the 
TOGA-COARE region from the SSM/I radiometer brightness temperatures (Tbs) 
have only moderate agreement with rain observations made by the ship-borne 
radars. In addition, these results indicate that mean rain rate deduced from these 
algorithms, representing averages for all the rain events observed by SSM/I over a 
time period of about 30 days (time span of a TOGA-COARE ship cruise) and over an 


1 Radiometer of the Defense Meteorological Satellite Program (DMSP) that has 19, 37 and 85 GHz 
channels in dual polarization and a 22 GHz channel in vertical polarization. This radiometer observes 
the earth's surface and atmosphere in a conical scan zvith an incidence angle of ~ 50 ' J . (for more details 
see Hollinger et al., 1985 ) 

2 Tropica! Ocean Global Atmosphere - Coupled Ocean Atmosphere Response Experiment 
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area of about 300x300 km 2 (the TOGA-COARE radar scan area) can differ 
systematically with respect to radar data. Because of such systematic errors, there can 
be an underestimation of 50 % in one month and an overestimation of 50 % in 
another month. Based on these findings, Prabhakara et al. (1998) conclude that the 
measurements made by the microwave radiometer do not contain sufficient 
information to retrieve good quality rain rates on a footprint scale. This conclusion 
is consistent with that of an earlier study by Schols et al. (1995). 

Rain retrieval on a footprint scale with microwave radiometers over the land 
is complicated further by terrain, vegetation, and surface wetness. The 
NASA/GSCAT method of Adler et al. (1994), and the NOAA/SRL algorithm based 
on the work of Grody (1991) and Ferraro et al. (1995), are designed to estimate rain 
rate on a footprint scale over land. These methods are rooted in the concepts 
provided by theoretical models (see for eg. Wu and Weinman, 1984; Kummerow et 
al., 1989 and Smith and Mugnai, 1992). The GSCAT method works primarily on the 
premise that rain rate increases as 85 GHz brightness temperature decreases, due to 
scattering by ice particles. The NOAA/SRL method is analogous, except it 
incorporates more of the data from the lower frequency channels to screen surface 
contamination. 

In a recent study, Conner and Petty (1998) find that the NASA/GSCAT and 
NOAA/SRL methods have limited success in retrieving rain information on a 
footprint scale over land. Furthermore, they show that these methods yield the 
highest Heidke skill score for detecting rain when the threshold on the rain rate is 
small (~1 mmhr' 1 ). That skill score decreases as the threshold rain rate increases. 
They find that the correlation coefficient between rain rate deduced from radar and 
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that retrieved from these methods on a footprint scale using the SSM/I data for a 
period of nine months over geographic areas of about 3°x3° is about 0.55. 

The microwave radiometer observations represent a vertical integral of the 
information about hydrometeors contained in a column of the atmosphere. Thus, 
the information about rain contained in the radiometer data is coupled in a non- 
linear way to that of cloud liquid water, and frozen and melting ice particles of 
different sizes, shapes, and densities (Schols et al., 1995 and Meneghini, 1996). 
Multi-channel, dual-polarization measurements of the microwave radiometer have 
a significant amount of redundant information (Prabhakara et al., 1992 and 
McGaughey and Zipser, 1996). Hence, the meager amount of independent 
information contained in the radiometer data cannot resolve all hydrometeors and 
their vertical distributions. For this reason, in the rain retrieval theoretical models 
the vertical profiles of hydrometeors above the cloud base are assumed. In 
particular these models assume that the density and particle size of frozen 
hydometeors are correlated with the rain rate produced by the clouds. A recent 
study of Houze (1997), based on aircraft Doppler radar observations of hydrometeors 
and winds in the TOGA-COARE experiment, indicates that this assumption is not 
true. This weakness of the algorithms resulting from the limited information 
content of the microwave radiometer leads to a poor quality of the rain rates 
deduced on a footprint scale. 

In this study, we assume that the amount of ice of all types in clouds is at a 
minimum and rain rate is near zero in a given pixel when the 85 GHz scattering 
depression indicated by the SSM/I data has a value near a threshold. Utilizing this 
threshold on the 85 GHz scattering, we can determine the fractional area enclosed by 
rain in a mesoscale region. Then, we have developed a rain rate retrieval technique 
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applicable to a mesoscale region that utilizes exclusively this fractional rain area. 
Thus, this method does not use the complete quantitative information of the 
85 GHz scattering depression on a footprint scale. 


We may remark that radar with its vertical resolving capability offers distinct 
advantages over the radiometer. Furthermore, radar back scatter measurements 
have a d 6 dependence on the rain drop diameter (Battan, 1973), while the 
radiometer observations tend to have about d 3 dependence (Olson, 1996). For these 
reasons, we contend radar data can provide good measurements of rain when 
suitably calibrated with surface rain gauges. However, the radar in the TRMM 3 
mission has less than a third of the scanning capability of the microwave radiometer 
(Simpson et al., 1996). This much larger scanning capability of the radiometer could 
be exploited to give more temporal and spatial coverage of rain information. 
However, we need a mesoscale rain retrieval technique based on the radiometer 
data that can be tuned with a limited amount of ground-truth radar observations. 
Such a method was developed for the SSM/I data on ocean (Prabhakara et al., 1998). 
The aim of this study is to develop this type of method for the land. 

2. Information about rain contained in the SSM/I data 

In order to understand in detail the relationship between the microwave 
radiometer Tb and observed rain rate, we make comparisons of the spatial 
distributions of brightness temperature in the 85 GHz horizontal polarization 
channel Tg 5 h and observations of rain rate for two mesoscale rain events. We are 
presenting one such comparison over ocean, and another over land to illustrate the 


3 Tropical Rainfall Measuring Mission (see Simpson et al., 1996).. 
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similarities in the microwave radiometer information over these regions. Because 
of the relatively small footprint size (~ 15 km) and strong extinction, the brightness 
temperatures in the 85 GHz channel are emphasized in these comparisons. 

Over the ocean, a rain event is chosen from the TOGA-COARE experiment, 
where observations of rain rate are made by ship-borne radars. These radars were 
calibrated with rain gauges on buoys and ships (Short et al., 1997). The SSM/1 data 
are matched with the radar data to within the time taken to complete a radar 
volume scan, i.e., 10 minutes. Furthermore, spatial matching between the two sets 
of data is within a few kilometers. 

As a land case, we have chosen a Mesoscale Convective System (MCS) event 
that occurred over the Northeastern United States in a 2°x3° grid box centered 
close to 79° W and40°N (the state of Pennsylvania). The distribution of observed 
rain rate on land in our study is obtained from 15-minute rain gauge data. About 
83 % of these rain gauges over the United States are Fisher-Porter (FP) rain gauges. 
These rain gauges have a poor sensitivity and can indicate only 15-minute 
accumulation of rain in increments of 0.1 in. In the remaining gauges, 15-minute 
accumulation of rain is measured in increments of 0.01 in. For uniformitv in this 
study, these rain rates are converted to mmhr 1 , and then rounded to the nearest 
unit. We find observations from about 40 rain gauges in a grid box of 2°x3° can 
lead to a representative average in that box. Such a density of gauges is present in 
some of the Northeastern United States. This density decreases in the Midwest, and 
is poor (~ 5 gauges per 2° x 3° grid box) in the Western states. For this reason, this 
study of precipitation on land is limited to a small section of the Northeastern 
United States, where there is a dense network of rain gauges. 


D 



The observations made by surface radars on land were not considered in our 
study. As pointed out by Conner and Petty (1998), these radar data have significant 
problems such as ground clutter, beam spreading, blockage by structures and 
topography, bright band, and anomalous propagation. They infer that because of 
several problems, surface radars can underestimate rain rate by as much as a factor 
of five. Hence, the data from rain gauges are commonly used to adjust the radar 
measurements (see for e.g. Oki et al., 1997). 

In Figure la, we show a map of Tssh over the TOGA-COARE region for a rain 
event that occurred at 1031 UTC on 24 Dec. '92 (Day 359). A rain rate map deduced 
from the ship-borne radars corresponding to this event is shown in Figure lb. Some 
regions with intense rain rate (> 4 mmhr 1 ) are identified in Figure lb with light 
shading; while other such intense rain regions are identified with cross hatching. 
This differentiation is done on the basis of the information given in the map of Tssh 
presented in Figure la. Generally, the light-shaded areas correspond to regions 
where Tssh is relatively warm (~ 260 K); while the cross-hatched areas correspond to 
regions where Tssh is cold (< 200 K). 

In Figure lc, a scatter plot of radar rain rate versus radiometer Tssh is shown 
to summarize the characteristics of the rain rates described above. We observe from 
this figure one maximum in the rain rate when Tssh is approximately 260 K, and 
another maximum when Tssh is around 200 K. These intense rain rates associated 
with relatively warm (cold) Tg 5 h are referred to as "warm" ("cold") convective rain. 
Only a few footprints show intense rain, while a large number near 260 K or 200 K 
show rain rates that are weak (approaching near zero value). These weak rain rates 
are broadly referred to as stratiform rain. From Figure lc, we note that weak rain 
rates can be present for all possible values of Tssh- 
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Figure 1: Maps of a) SSM/I Tg5h and b) the coincident radar rain rate data in the TOGA-COARE region for 24 Dec. '92 at 1031 
UTC. In lb, regions with intense rain rate (> 4 mmhrl) are identified with light shading and cross hatching. This ditferentiation 
is done on the basis of the information given in the map of Ts5h presented in la. Light-shaded areas correspond to regions where 
Tg5h is relatively warm (~ 260 K); while the cross-hatched areas correspond to regions where Tg5h is cold (< 200 K). c) Scatter 
plot of coincident 85 GHz SSM/I radiometer brightness temperature (K) measurements and radar rain observations (mmhr^) over 
the TOGA-COARE region, d) Map of SSM/I T37h in the TOGA-COARE region for 24 Dec. '92 at 1031 UTC. e) Scatter plot of 
coincident 37 GHz SSM/I radiometer brightness temperature (K) measurements and radar rain observations (mmhr*) over the 
TOGA-COARE region. 





In Fig Id a map of the 37 GHz SSM/I data corresponding to the above TOGA- 
COARE rain event is shown to illustrate the redundancy in the information 
conveyed by the 85 and 37 GHz channels. In several large areas (> 50x50 km 2 ), where 
the 85 GHz channel shows a significant depression, the 37 GHz channel shows a 
noticeable warming. Furthermore, the scatter plot shown in Figure le of T 37 h 
versus rain rate for this event indicates that near zero rain rate (averaged over the 
37 GHz footprint) can be present at all possible values of T 37 h- This property of T 37 h 
is similar to that of T§5h shown in Figure lc. Thus, we find from these observations 
that estimation of rain rate for given values of Tg 5 h and/or T 37 h can be totally non- 
representative of the ground-truth rain rate. This is also true of the 19 and 22 GHz 
channels. 

In Figures 2a, b, c, and d, we are presenting information pertaining to an MCS 
event that took place over the Northeastern United States in a 2°x3° grid box 
centered close to 79° VV and40°N. This event happened at 619 EST on 22 Aug. '87. 
In Figure 2a, we present an analysis of Tgsh data that is similar to that in Figure la. 
In Figure 2b, we present the distribution of rain rate on land for this event using the 
available data from rain gauges. We may emphasize that the heavily shaded area of 
Figure 2a, where strong 85 GHz scattering is observed (Tssh < 200 K), is at least 
50x100 km 2 . Although this is the region with the heaviest scattering, the maximum 
rain intensity of 30 mmhr 1 is observed 50 km away in a region where Tssh is about 
240 K. Furthermore, rain intensity varies somewhat randomly from 0 to 10 mmhr 1 
over the area encompassed by the 230 K contour of Tg 5 h- We infer from these two 
figures that the rain rate on land does not necessarily increase as Tggh decreases. 
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On land, when the surface gets wet, its emissivity decreases and hence the 
37 GHz brightness temperature (T 37 h) is diminished as compared to that on dry land 
(Prabhakara et al, 1995). This reduction in T 37 h can be noticed because the radiation 
at this frequency is weakly absorbed by water vapor and oxygen in the atmosphere. 
In the 85 GHz channel of SSM/I, this absorption due to water vapor and oxygen in 
the atmosphere is relatively strong. For this reason, the effect of the surface wetness 
on Tgsh is much less pronounced. In Figure 2c, we show a map of the difference 
signal, T 37 h - Tgsh- This difference is used to minimize impact of the surface 
wetness effect on Tssh- Essentially all of the salient features in Figure 2a are 
represented here. We may note that because of scattering, the brightness 
temperature in the 85 GHz channel is less than that of the 37 GHz. Thus, the 
difference signal, T 37 h - Tggh, serves as a scattering index on land (see also Barrett et 
al., 1988), and is used in the subsequent rain analysis. 

In order to show information about rain in the long wavelength channels at 
19 and 37 GHz, we show in Figure 2d a map of the difference Tj 9 h - T 37 h for this rain 
event. Examination of Figures 2a, c, and d reveal that there is substantial 
redundancy in the information of the 85, 37, and 19 GHz channels. However, 
information given by the high frequency channels (T 37 h and Tgsh) is preferred 
because of better sensitivity and spatial resolution. 

Since the rain gauges have poor resolution, and are not dense in spatial 
coverage, we cannot show a scatter plot of rain rate versus Tg5h for one event, as 
was done in Figure lc. For this reason, we have taken the 15-minute rain gauge data 
over a large area of the Eastern United States (25-45 N and 75-84 W) for the month of 
Aug. '87 and compared it with the SSM/I Ts5h measurements. We may note that 
these gauge and radiometer data are matched in space and time to within 4 km and 
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15 minutes, respectively. In Figure 3a, we show a scatter plot of these data. A 
similar figure for the Midwest United States (25-45 N and 88-97 W) is presented in 
Figure 3b. We notice from these plots that the rain rate is not related to Tssh in a 
simple fashion, which was found to be the case in Figure lc over ocean. In 
Figure 3a, a prominent maximum near 230 K is apparent; while in Figure 3b, the 
"warm" and "cold" convective rain maxima are clearly seen. From these analyses 
shown Figures la-e, 2a-d, and 3a-b, one can see that the signal due to rain in the 
microwave radiometer data on land and ocean is broadly similar. 

We may remark that the SSM/I data in the 19, 37, and 85 GHz channels are 
measured in dual polarization. We have examined the polarization information in 
all these channels, together with the rain observations over land and ocean. We 
find addition of this polarization information does not necessarily help us to 
improve quantitative estimation of rain rate. 

From the above discussion, we find that the information contained in the 
multi-spectral, dual-polarization measurements of SSM/I is highly redundant. 
When there is rain, other hydrometeors - cloud liquid water and frozen and melting 
ice particles of different sizes, shapes, and densities - can also be present in a vertical 
column of the atmosphere. The radiometric contamination introduced by these 
other hydrometeors is not necessarily insignificant. Thus, with the limited 
information from the microwave radiometer, theoretical rain retrievals (see for eg. 
Wu and Weinman, 1984; Kummerow et al., 1989 and Smith and Mugnai, 1992) 
encounter insurmountable problems. We may note that this is not solely because of 
beam-filling and surface emissivity problems. 
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Figure 3: Scatter plot of coincident SSM/I Tg5h vs. rain gauge observations (mrnhr'l) over land in a) the Eastern United States 
(25-45 N and 75-84 W) and b) the Midwest United States (25-45 N and 88-97 W). 



Because of these problems, in this study we have taken the path of mesoscale- 
average rain rate estimation using the microwave radiometer data. The basis for 
this approach is that a threshold on the 85 GHz scattering depression can be used to 
detect the presence of rain. This allows us to determine the rain area in a mesoscale 
region that exceeds this threshold. Then, this rain area can be related to mesoscale- 
average rain rate. By doing so, we ignore the details of the nature of rain on a small 
scale; but we attempt to capture the intensity of rain on a macroscopic scale. 

3. Area-average rain retrieval method applicable to land 

In this study, we derive the fractional rain area for three individual grid boxes 

of 2°x3° over a section of the Northeastern United States (39.5-41.5 N and 

81-75 W) with the aid of T37h - Tssh/ which is denoted as 6T37-85- These three grid 

boxes exclude the area of the Great Lakes. The fractional rain area in a mesoscale 

grid box of 2° x 3° is indicated with the symbol fR, where fR is given by 

, n 

r- n' (1) 

In this equation, n is the number of 85 GHz rain footprints in a 2° x 3° grid box, and 
N is the total number of 85 GHz footprints in that grid box. We have used values of 
1, 2, 3, 4, and 8 K for 5T37.S5 as thresholds to detect the presence of rain. Of these 
thresholds, the best linear correlation between fR and average rain rate is found 
when 5T37-85 is 3 K. The 3 K threshold on 8T37.85 is used throughout the remainder 
of this study to determine fR. 

We have compiled nine months of SSM/I measurements (Mar. '88 - 
Nov. '88) in the three grid boxes. In this compilation, only those passes that have 
grid boxes that are fully covered by the satellite data are retained. From these data, 
we have deduced fR for each grid box and the corresponding mean of ST37-85 over 
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the rain area fR , i.e., ( 8X37-85 )■ In a given grid box for each satellite pass, the 
ground-truth, grid-box-average rain rate is obtained by averaging the 15 -minute rain 
gauge data. In order to ensure a robust estimate of this average, we require that 
there should be at least 40 gauges in a given grid box. This average rain rate is 
denoted by R>4o- This R>40 is compared with fR and ST 37-85 with the data 
presented in Table 1 for one grid box ( 39 . 5 - 41.5 N and 78-81 W) for the full nine 
months. 

In Figure 4 a, a scatter plot of R>40 versus fR is presented for the nine months 
of data given in Table 1. The correlation coefficient between these variables is 
shown in the figure. Similarly, in Figure 4 b, a scatter plot of R>40 versus 6T37_g5 
and the correlation coefficient between these two variables is presented. From these 
figures, we note that fR explains a significantly larger fraction of the variance in 
R>40 than 5T37_85. 

We have computed grid-box-mean, i.e. mesoscale-average, rain rate using a 
retrieval formula that is similar to that presented in Prabhakara et al. ( 1998 ). The 
form of that equation is given below: 

R X =ex Pkx X ]- 1 (2) 

The retrieved rain rate Rx given by the above equation depends on the argument X 

and the scaling constant vyx shown in the exponent. Substituting the argument X 
with fR and ST 37-85, we have estimated mesoscale-average rain rates Rf and 

R 37 - 85 / respectively. In Figures 5 a and 5 b, R>40 is compared against Rf^and R37_s5, 
respectively. We find that Rf^ explains 66 % of the variance contained in R>40, 

while R37-85 explains only 21 %. Although we have not presented for this case 
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Figure 4: a) Relationship between the fractional rain area fR and the corresponding grid-box-average rain rate R > 4 o deduced from 
15-minute rain gauge observations, b) Relationship between the scattering index (T37h-T§5h) and the corresponding grid-box- 
average rain rate R> 40 deduced from 15-minute rain gauge observations. 
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Figure 5: Scatter plot of R> 4 q versus grid-box-average rain rate retrieved using the equations a) R tR =exp[0.036*f R j-l and 

b) ^ 37.85 = exp| 0.003 • f R • 5T 37 - 35 ] -1. These plots represent data for nine months of 1988 (Mar. '88 - Nov. ' 88 ) in a grid box (39.5- 
41.5 N and 78-81 W) in the Northeastern United States. 




results from substituting X with simple linear and non-linear combinations of the 
two basic variables fR and 6X37-85, we find that the retrieved rain rates from such 
combinations explain significantly less than 66 % of the variance contained in R>40- 

The NASA/GSCAT and NOAA/SRL rain retrieval algorithms, applicable to 
the radiometer footprint scale, are presented in detail by Conner and Petty (1998). 
These algorithms have built in thresholds on the 85 GHz brightness temperature, as 
well as some screening procedures to eliminate questionable surface emissivity 
conditions. We have implemented these algorithms in this study and computed 
rain rates on a pixel scale for each satellite pass in each grid box. In order to compare 
these results with R>4o, we have averaged these pixel-scale rain rate estimates in 
each grid box for each satellite pass. These grid-box average rain rates for the 
NASA/GSCAT and NOAA/SRL rain retrieval algorithms are denoted as Rgscat 
and Rsrl, respectively. We may remark that since both Rgscat ar *d Rsrl are 
based on thresholds on the 85 GHz brightness temperature, these grid-box average 
rain rates reflect the joint effect of the rain area in the grid box as well as the 85 GHz 
scattering index. In Figures 6a and 6b, we have compared these rain rates with R>40- 
The NASA/GSCAT and NOAA/SRL rain retrieval algorithms can explain 
respectively about 18 % and 14 % of the variance contained in nine months of the 
ground-truth rain rate data over this grid box in the Northeastern United States. 

We have retrieved grid-box-average rain rates for the other two grid boxes 
over the Northeastern United States for the same nine-month time period. These 
are the only two other 2° x 3° grid boxes in the United States where the number of 
15-minute rain gauges exceeds 40 in each grid box. In Table 2 combining the data for 
all three grid boxes, we compare the relative merits of the rain rate retrieved from 
Equation 2 when X is substituted by fR and the product fR *6X37-85 • The rain rates 
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Gauge vs. NASA/GSCAT Algorithm 



Grid-Box-Average Rain Rate: NASA/GSCAT Algorithm (mmhr' 1 ) 
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Grid-Box-Average Rain Rate: NOAA/SRL Algorithm (mmhr 1 ) 


Figure 6: Scatter plot of R> 4 q versus mesoscale-average rain rate retrieved using a) NASA/GSCAT and b) NOAA/SRL 

algorithms. These plots represent data for nine months of 1988 (Mar. ’88 - Nov. ’88) in a grid box (39.5-41,5 N and 78-81 Wi in the 
Northeastern United States. 



retrieved from the NASA/GSCAT and NOAA/SRL algorithms are also included in 
this table for comparison. From this comparison, we conclude that the grid-box- 
average rain rate retrieved from fR , shown in the first row of the table, can explain 
about 64 % of the variance contained in the ground-truth rain observations in all 
three grid boxes. Performance of this algorithm is better than that of the others, 
except in one grid box. 

The reason why the algorithm based on fR performs better than the 
NASA/GSCAT and NOAA/SRL algorithms is explained as follows. In the 
fractional rain area algorithm, only fR gives the grid-box-average rain rate. On the 
other hand, in the other algorithms both the rain area and the scattering 
information enter into the retrievals. As long as fR and the scattering information 
are correlated, we expect good retrievals. On the contrary, when this correlation is 
poor, the retrievals are degraded. The correlation between fR and the scattering 
information in a given period and grid box deteriorates when there are some events 
in that period where fR is relatively large and the scattering information is weak. 
The poor results shown in rows 3 and 4 of Table 2 for grid boxes 2 and 3, 
representing the performance of the rain retrievals of NASA/GSCAT and 
NOAA/SRL, can be explained with the above reasoning. 

When the data from all three grid boxes are put together, we find that Rf^can 

explain 50 % more of the variance contained in R> 40 , as compared to RgscAT and 
Rsrl • All three algorithms have inherent in them the concept of the rain area. To 
a first-order approximation, one may expect that the average rain rate in a grid box 
increases as the rain area increases. However, because the NASA/GSCAT and 
NOAA/SRL algorithms use additionally the information of the 85 GHz scattering 
index, these rain retrievals deteriorate. This result can be inferred from the non- 
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monotonic relationship between the ground-truth rain rate and Tssh shown in 
Figures lc, 3a and 3b on a footprint scale. Because of this reason, even when the 
ground-truth rain rate and Tgsh are averaged over a mesoscale area, they relate only 
poorly with one another. 

4. Procedure to retrieve rainfall on land in the TRMM mission 

The method described earlier gives a radar-quality, monthly-mean rain rate 
in a mesoscale grid box. A brief summary of the procedure to apply this rain 
retrieval method to TRMM radiometer and radar data is presented in the following: 

1) From the 760 km wide swath of TRMM Microwave Imager (TMI) radiometer 

data, extract a 220 km swath of data centered on the satellite orbital track. This 
data overlaps exactly with the TRMM radar observations. Bin these TMI data 
into a fixed geographic grid of 2° x3°. For each pass of the TRMM satellite 
over a given grid box, the collocated TMI and radar data will generally fill a 
variable portion of that grid box. This variable portion of a grid box is 
considered to be a mesoscale area. For purpose of TRMM, we will accept for 
the rain retrievals only those events that fill at least 50 % of a given grid box 
with collocated TMI and Precipitation Radar data. 

2) Within a given month, for each grid box in the land regions from 35 S to 35 N 

covered by the TRMM satellite, do the following: 

a) Determine the number of rain events. Then, utilizing the TMI data, 
estimate the fractional rain area, f^, for each of these rain events. Also, 
from the radar data, determine the average rain rate Ra for each rain 
event. 
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the radiometer and the sum of Ra from radar for all the rain events in 
the month are equal. 


N , x N 

X( R f J. = 51 ( R A)j • 

i=i v k a i.i 

In the above summation, N m stands for the total number of rain events 

in that month. The above condition also means that the monthly- 
mean values Rf and Ra are equal. 

Procedures 1 and 2 will create a map of with a resolution of 2°x3° over 

R 

the land areas from 35°N to 35°S for a given month. Then, using all of the 760 km 
wide swath of TMI radiometer data, we can derive area-average rain rates in each 

2° x 3° grid box for each rain event detected by the TMI. In this manner, in a month 
each grid box is sampled by the TMI about 3.5 times more often than by the radar. 
From these gridded data, the 5°x5° monthly-mean rainfall as required by the 
TRMM mission (Simpson et al., 1996) can be determined by suitable averaging. 

5. Discussion and conclusions 

In this study, we indicate that on a footprint scale the limited amount of 
information contained in the multi-channel, dual-polarization SSM/I microwave 
radiometer measurements over land and ocean is not sufficient to resolve different 
types of hydrometeors or rain intensity in a satisfactory fashion. This is primarily 
because the information conveyed by these channels is highly redundant and is not 
capable of discriminating different kinds of hydrometeors and their vertical profiles. 
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The radiometer data on the other hand can identify reasonably well the 
presence of rain for rain rates above a relatively small threshold. This allows us to 
determine the rain area in a mesoscale region that exceeds this threshold. We find 
that as this rain area increases the mesoscale-average rain rate increases. This result 
is simple and is consistent with the findings of earlier investigators (Doneaud et al., 
1984, and Lopez et al., 1989). Area-average rain rate in three 2° x3° grid boxes over 
the Northeastern United States (39.5-43.5 N and 75-81 W) estimated with this 
retrieval method for nine months has a correlation of 0.79 with rain rate derived 
from 15-minute rain gauges. We find that adding the spectral information from the 
radiometer to the retrieval algorithm causes the results to deteriorate. This result is 
borne out by the NASA/GSCAT and NOAA/SRL algorithms, where inherently the 
rain area and spectral information are combined. This result can be inferred from 
the non-monotonic relationship between the ground-truth rain rate and Tg5h 
shown in Figures lc, 3a and 3b. 

A similar analysis (Prabhakara et al., 1998) presented for the TOGA-COARE 
region showed that the spectral information helped only marginally in explaining 
the variance of the ground-truth radar rain rates. However, on the land, probably 
because of more vigorous convective activity, we find that the spectral information 
degrades the rain retrievals. 


We find that the area-average rain retrieval method requires ground truth 
that has good sensitivity and areal coverage. This requirement will be met when 
the TRMM radar rain rate observations are available. Additionallv, the field of the 
view of the TRMM radiometer is about one third that of the SSM/I, which should 
enable us to define the fractional rain area much more accurately. We believe that 
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these benefits should give improved results when the algorithm is applied to the 
TRMM data. 


In the TRMM mission, the radar observations are limited to a narrow swath 
of 220 km, while the microwave radiometer observations extend over a swath of 
760 km. Using the coincident TRMM radiometer and radar data, we can calibrate 
our area-average rain rate retrieval method. With this calibrated method, we can 
extend the temporal and spacial coverage of the radar observations, and thus get an 
improved rainfall climatology. 
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Figure Captions 


Figure 1: Maps of a) SSM/I Ts 5 h and b) the coincident radar rain rate data in the 
TOGA-COARE region for 24 Dec. '92 at 1031 UTC. In lb, regions with intense rain 
rate (> 4 mmhr 1 ) are identified with light shading and cross hatching. This 
differentiation is done on the basis of the information given in the map of Tssh 
presented in la. Light-shaded areas correspond to regions where Tssh is relatively 
warm (~ 260 K); while the cross-hatched areas correspond to regions where Tgsh is 
cold (< 200 K). c) Scatter plot of coincident 85 GHz SSM/I radiometer brightness 
temperature (K) measurements and radar rain observations (mmhr 1 ) over the 
TOGA-COARE region, d) Map of SSM/I T 37 h in the TOGA-COARE region for 24 
Dec. ’92 at 1031 UTC. e) Scatter plot of coincident 37 GHz SSM/I radiometer 
brightness temperature (K) measurements and radar rain observations (mmhr -1 ) 
over the TOGA-COARE region. 


Figure 2: Maps of a) SSM/I Tssh/ b) rain rate from 15-minute rain gauge, c) SSM/I 

T37h-T85h/ and d) SSM/I Ti 9 h-T 37 h over a 2°x3° grid box in the Northeastern 
United States on 22 August '87 at 1119 UTC (619 EST). In 2a, the prominent regions 
of strong scattering depression, i.e. Tssh < 200 K, are denoted by heavy shading. The 
regions where Tssh is between 200 K and 230 K are lightly shaded. 

Figure 3: Scatter plot of coincident SSM/I Tgsh vs. rain gauge observations 
(mmhr -1 ) over land in a) the Eastern United States (25-45 N and 75-84 YV) and 
b) the Midwest United States (25-45 N and 88-97 VV). 



Figure 4: a) Relationship between the fractional rain area f^ and the corresponding 
grid-box-average rain rate R >40 deduced from 15-minute rain gauge observations, 
b) Relationship between the scattering index (T37h-Ts5h) and the corresponding grid- 
box-average rain rate R >40 deduced from 15-minute rain gauge observations. 


Figure 5: Scatter plot of R >40 versus grid-box-average rain rate retrieved using the 
equations a) Rf R = exp[0.036 • f R ] - 1 and b) R 3 „ 8 _ =exp[o.003 * f R •5T 37 _ 85 ]-1. These 


plots represent data for nine months of 1988 (Mar. '88 - Nov. ’88) in a grid box 
(39.5-41.5 N and 78-81 YV) in the Northeastern United States. 


Figure 6: Scatter plot of R> 4 o versus mesoscale-average rain rate retrieved using 

a) NASA/GSCAT and b) NOAA/SRL algorithms. These plots represent data for 
nine months of 1988 (Mar. ’88 - Nov. '88) in a grid box (39.5-41.5 N and 78-81 W) in 
the Northeastern United States. 



Table 1 . The data in a 2°x3° grid box within the Northeastern United States 
(39.5-41.5 N and 78-81 W) for the nine month period Mar. '88 to Nov. '88 used for 
area-average rain retrieval. The data is as follows: 

Jday/Mon - Julian day and month 

Pass - Satellite pass relative to 0000 UTC of the 1st day of a given month 
n - Number of 85 GHz rain pixels 
N - Total number of 85 GHz pixels 
fR - n/N Fractional rain area 


T37-85av - 5T37-85 averaged over fR, i.e., 5T 37.35 
#G - Number of rain gauges 


RG - Area-average rain rate deduced from 15-minute rain gauges, i.e., R>40 
R1 - Area-average rain rate retrieved using Rf R = exp[0.036* f R ]-l 


R2 - Area-average rain rate retrieved using R 37 g . =exp|o.003 • f R • 5T 


37-85 


-1 


R3 - Area-average rain rate retrieved using NASA/GSCAT algorithm 
R4 - Area-average rain rate retrieved using NOAA/SRL algorithm 
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Table 2. Correlation coefficients between retrieved grid-box-average rain rate and 
gauge rain rate for three grid boxes in the Northeastern United States for nine 
months of 1988 (Mar. ’88 - Nov. ’88). 


Rain Retrieval Method 


Grid Box 1 Grid Box 2 

(39.5-41.5 N, (39.5-41.5 N, 

78-81 W) 75-78 W) 


R f R =ex P[0 036 • f r ] — 1 0.81 


0.79 


R 


37-85 


0.64 


0.85 


expjo.003 • fR • ST 37 _ 85 


NOAA/GSCAT 


0.43 


0.78 


Grid Box 3 
(41.5-43.5 N, 
75-78 W) 

0.81 


0.73 


0.67 


NOAA/SRL 


0.37 


0.89 


0.67 
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